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Abstract 

We report a precise measurement of the Df meson lifetime. The data were taken 
by the SELEX experiment (E781) spectrometer using 600 GeV/c E - , ir~ and p 
beams. The measurement has been done using 918 reconstructed Djr. The lifetime 
of the Df is measured to be 472.5 ±17.2 ±6.6 fs, using K*(892)°K ± and 0vr ± decay 
modes. The lifetime ratio of Df to D° is 1.145 ± 0.049. 
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1 Introduction 



Precise measurements of the lifetimes of charm meson weak decays are im- 
portant for understanding QCD in both perturbative and nonperturbative 
regimes. For mesons a joint expansion in Heavy Quark Effective Theory and 
perturbative QCD parameters treated through the third order in the heavy 
quark mass shows a term including non-spectator W- annihilation as well as 
Pauli interference. The resulting non-leptonic decay rate differences between 
W-exchange in D° and W-annihilation in Df produce lifetime differences of 
order 10-20 % (0). 

The Df lifetime (§ was dominated by the measurements from E687 Collab- 
oration (0.475 ± 0.020 ± 0.007 ps) (||). Recently new precision measurements 
of the Df lifetime have been made by the E791 Collaboration (0.518 ± 0.014 
± 0.007 ps) (D and the CLEO Collaboration (486.3 ± 15.0 ±f;? fs) (§). Both 
groups have taken advantage of improved precision in the D° lifetime measure- 
ment to report new results for the DftoD lifetime ratio of 1.254 ± 0.041 (^) 
and 1.190 ±0.042 (||). Their average is 7.4 a from unity, emphasizing the large 
difference in W contributions to Df and D° decays. 

In this letter we report the results of a new measurement of the Df lifetime 
based on data from the hadroproduction experiment SELEX (E781) at Fer- 
milab. The measurement is based on about 1000 fully reconstructed decays 
into KKtt from a sample of 15.3 xlO 9 hadronic triggers. 

The SELEX detector at Fermilab is a 3-stage magnetic spectrometer. The 
negatively charged 600 GeV/c beam contains nearly equal fractions of £ and 
7i. The positive beam contains 92% protons. Beam particles are identified 
by a Transition Radiation detector. The spectrometer was designed to study 
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charm production in the forward hemisphere with good mass and decay vertex 
resolution for charm momenta in the range 100-500 GeV/c. Five interaction 
targets (2 Cu and 3 C) had a total target thickess of 4.2% X int for protons. 
The targets are spaced by 1.5 cm. Downstream of the targets are 20 silicon 
planes with a strip pitch of 20-25 /zm oriented in X, Y, U and V views. The 
scattered-particle spectrometers have momentum cutoffs of 2.5 GeV/c and 
15 GeV/c respectively. A Ring-Imaging Cerenkov detector (RICH) (|6p, filled 
with Neon at room temperature and pressure, provides single track ring radius 
resolution of 1.4% and 2a K/ir separation up to about 165 GeV/c. A layout 
of the spectrometer can be found elsewhere (0). 



2 Data set and charm selection 

The charm trigger is very loose. It requires a valid beam track, at least 4 
charged secondaries in the forward 150 mrad cone, and two hodoscope hits 
after the second bending magnet from tracks of charge opposite to that of the 
beam. We triggered on about 1/3 of all inelastic interactions. A computational 
filter linked PWC tracks having momenta > 15 GeV/c to hits in the vertex 
silicon and made a full reconstruction of primary and secondary vertices in 
the event. Events consistent with only a primary vertex are not saved. About 
1/8 of all triggers are written to tape, for a final sample of about 10 9 events. 

In the full analysis the vertex reconstruction was repeated with tracks of all 
momenta. Again, only events inconsistent with having a single primary vertex 
were considered. The RICH detector identified charged tracks above 25 GeV/c. 
Results reported here come from a preliminary reconstruction through the 
data, using a production code optimized for speed, not ultimate efficiency 
The simulated reconstruction efficiency of any charmed state is constant at 
about 40% for x F > 0.3 where > 60% of SELEX events lie. 

To separate the signal from the noncharm background we require that: (i) the 
spatial separation L between the reconstructed production and decay vertices 
exceeds 8 times the combined error o^, (ii) each decay track, extrapolated to 
the primary vertex z position, must miss by a transverse distance length t > 
2.5 times its error <r t , (iii) the secondary vertex must lie outside any target by 
at least 0.05 cm and (iv) decays must occur within a fiducial region. 

There are 918 ± 53 events Df candidates, each having two RICH-identified 
kaons and a pion, for which no particle identification is required. We divide 
them into three decay channels: K*(892)°K ± , cf)^ and other KK7T. The res- 
onant mass window for the ir*(892)° (0) was 892 ± 70 MeV/c 2 (1020 ± 10 
MeV/c 2 ). 
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n/K misidentification causes a reflection of D ± under the Df peak. We limit 
the maximum kaon momentum to 160 GeV/c to reduce misidentification in the 
RICH. To evaluate the shape of this background we use the D ± — > K^ir ir 
sample that passes all the cuts listed above and lies within ± 15 MeV/c 2 of the 
D + mass. We formed the invariant mass distribution of these events when one 
pion is interpreted as a kaon. At most one of the two possible reflections per 
event falls into the Df mass window. The reflected mass distribution was fit 
by a polynomial function rising at 1925 MeV/c 2 and decreasing to zero at large 
invariant mass. Dividing this distribution by the number of D ± events gives 
us the contribution per mass bin for each misidentified D^ in the Df sample. 
We count the misidentified D in the Df sample by fitting the Df mass 
distribution within ± 20 MeV/c 2 interval around the Df mass with the sum 
of a Gaussian signal, a linear background shape estimated from the sidebands 
and the D shape with variable normalization. The resultant misidentified D^ 
contribution to the Df mass distribution is shown as the hatched areas in Fig. 
1(a), (b). The fit gives 52±7 and 12±3 misidentified D ± events in the K*K and 
(j)7i decay mode, respectively (the error quoted is statistical only). The RICH 
kaon identification is a very powerful tool for rejecting D^ contamination; <p 
decay kinematics further reduces particle identification confusion in the (pir 
channel. 

To estimate yields, we subtracted the sideband background and D ± contami- 
nation as evaluated above from the total number of events in the signal region. 
We find 430 ± 24 K*K and 330 ± 19 0tt events. A Gaussian fit to the com- 
bined data shown in Fig. 1(a), (b) gives a mass of 1969 ±9.7 MeV/c 2 and 
a D ± = 8.0 MeV/c 2 . Note that these fitted yields are not used as constraints 
in the lifetime fit, as discussed below. 



3 Lifetime evaluation using a maximum likelihood fit 

The average longitudinal error a z on the primary and secondary vertices is 270 
/im and 500 /im, which gives a combined error of 570 /im. In the Df sample, 
the average momentum is 215 GeV/c, corresponding to a time resolution of 18 
fs, about 4% of r D ± . Because bin-smearing effects are small, we used a binned 
maximum likelihood fitting technique to determine the Df lifetime. The fit was 
applied to a reduced proper time distribution, t* = M{L — L min )/pc where M 
is the known charm mass (|2]), p the reconstructed momentum, L the measured 
vertex separation and L min the minimum L for each event to pass all the 
imposed selection cuts. L min is determined event-by-event, along with the 
acceptance, by the procedure described below. We fitted all events with t* < 
1600 fs in the mass range 1.949 < M(KKn) < 1.989 GeV/c 2 , ± 2.5a from 
the Df central mass value. 
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To evaluate the mean lifetime we used a maximum likelihood method. The 
probability density was performed by the function : 



f(r Da , T B ck, a, /3; t*) = 

(1 - ^)N s e -— — + aN s B{t) + ^ e -*V'b± (1) 
e{t*)T D ± t d ± 

where 

g e -t*/T B ok l - R 
B(t*) = P - + — -£ (2) 

The function is the sum of a term for the Df exponential decay corrected 
by the acceptance function e(t*) plus a background function B(t*) consist- 
ing of a single exponential plus a constant to account for a flat background 
extending to large proper time. Its parameters were determined from the t* 
distribution from the Df sidebands. It also includes a term for the D ± expo- 
nential decay normalized to the number of misidentified events in the signal 
region. The t d ± lifetime used in the fit is 1051 ±13 fs (0). The mass range 
of the sideband background windows, 1.890 < M(KKn) < 1.930 GeV/c 2 and 
2.040 < M(KKn) < 2.080 GeV/c 2 was twice the signal mass window. We 
defined asymmetric sidebands to avoid the influence of — > K + K~ tt , 
and we excluded the £)*(2010) mass region. 

The four parameters are: r D ± (Df lifetime), TBck (background lifetime), a 
(background fraction in the signal region) and (3 (background splitting func- 
tion). Ns is the total number of events in the signal region after D contam- 
ination subtraction. 



The proper-time-dependent acceptance e(t*) is independent of spectrometer 
features after the first magnet, e.g., RICH efficiency and tracking efficiency 
These efficiencies affect only the overall number of events detected. The proper 
time distribution of these events depends crucially on vertex reconstruction. 
To evaluate e(t*) we reanalyze each observed Df event after moving it to a 
large set of different proper times t*. Only the longitudinal position of the 
charm decay point and the axial orientation of the 3-body decay vectors are 
changed (H) . A reanalyzed event is accepted if it passes the same cuts as those 
applied to the data. The minimum flight path after which an event is accepted 
defines L m i n for this event. This method is independent of details of the true 
Xf or transverse momentum distributions (§). 

Fig. 2 shows the overall fits to the data distributions as a function of reduced 
proper time for K*(892)°K ± and (pir^ decay modes. It also shows the accep- 
tance e(t*), which differs significantly from unity only after 4 lifetimes where 
statistics are limited. 
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Df 


r D ±(fs) 


Signal 


^(892)°^ 


472.3 ± 23.0 


430 ± 24 




473.0 ± 26.0 


330 ± 19 


Average 


472.5 ± 17.2 


760 ± 30 



Table 1 

Lifetime results and signal yields for the two Df modes analyzed. The last row 
is the weighted average of the two resonant channels K* (892)° and cftir^. The 
errors are statistical only. 



Source of uncertainty 


K*(892)°K ± 


07T ± 


Vertex reconstruction 


<1 


<1 


D contamination 


2 




Acceptance function 


3 


2 


Fit procedure 


5 


0.3 


Total systematic error 


6.2 


2.2 



Table 2 

Systematic error contributions in fs. 

Table 1 summarizes the lifetime results for the two modes analyzed: Df —> 
K*(892)°K ± ; and Df — ► (pn^. The uncertainties are statistical only, evaluated 
where — ln£ increases by 0.5. Combining these results for the two resonant 
modes, we measure an average lifetime t d ± = 472.5 ± 17.2 fs. 



4 Systematic errors 

The systematic uncertainties for the Df lifetime analysis are listed in Table 2 
and described below. We group them in the following categories: 

4-1 Primary and secondary vertex reconstruction 

Lifetime shifts due to reconstruction errors have been well studied in our D° 
and A c work, with an order of magnitude higher statistics (§; ||). Because of 
the high redundancy and good precision of the silicon vertex detector, vertex 
mismeasurement effects are small at all momenta. Proper time assignment 
depends on correct momentum determination. The SELEX momentum error 
is less than 0.5% in all cases. We assign a maximum systematic error from 
proper time mismeasurement of 1 fs. 
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4-2 Misidentification 



The effect of contamination under the Df peak was studied by changing 
the width of the exclusion window around the nominal D ± mass for the K/n 
interchange discussed above. Effects on the (pir mode are negligible. For the 
K*{892)°K ± mode this gives a systematic error of 2 fs. 

4-3 Acceptance function 

The technique to determine the acceptance correction dependence on proper 
time is discussed extensively in Ref . (H) . It has been verified with much larger 
statistics there. The maximum systematic error here is dominated by the 
K*(892)°K ± correction, 3 fs. For the (pir^ 1 mode it is less than 2 fs. 

4-4 Fit procedure 

The fit was performed by the maximum likelihood method using a background 
parametrized by an exponential function plus a constant. We varied the width 
of the sidebands and the t* Max independently. The sytematic error due to the 
fit procedure is 5 fs and less than 0.5 fs for K*(892)°K ± and ifiir^ decay 
modes respectively. That error is mainly dominated by the Df background 
parametrization. 

Combining in quadrature all the sources of systematic errors listed in Table 2 
we obtain a total systematic error of 6.2 fs (2.2 fs) for the K*(892)°K ± (07r ± ) 
mode. 



5 CONCLUSIONS 

We have made a new measurement of the Df lifetime in two independent 
resonant decay channels, K*(892)°K ± and (ftir^ 1 using a maximum likelihood 
fit. SELEX measures the Df lifetime to be 472.5 ± 17.2 ± 6.6 fs. Using the 
result reported in the PDG (§) t d o = 412.6 ± 2.8 fs we evaluate a ratio 
t d ±/t d o = 1.145 ± 0.049, 3a from unity. The D s lifetime reported in this 
Letter is comparable in precision with previous experiments (f|; |5|). Our result, 
combined with other world data, (|2p, lowers the overall Df lifetime somewhat. 
Nevertheless, it is clear that the lifetime ratio t d ± / r^o is significantly larger 
than unity for all the precision measurements. 
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Fig. 1. Invariant mass distributions for a) — > i^*(892)°i^ :l= , b) Z)^ — > (^Tr 1 * 1 . 
Hatched regions show the K + K~ tt^ background from misidentified -D 1 * 1 . The arrows 
show the Df signal region. 
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Fig. 2. Corrected reduced proper time distributions for events in the Dg window 
1969 ± 20 Mev/c 2 (full circles) and results from the maximum likelihood fit (solid 
curve for: a) D s ± -> K*(892)°K ± ; and b) D s ± -» 0^ ± . The dashed curve shows the 
fitted background, including the contribution. The dashed-dot curve describes 
the acceptance function, e(t*). 
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